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Motivation 


Paths to transition 1 : The rule of transient growth 


» The non-normality of the NS equations can 
give rise to transient energy amplification 
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9 Transient growth is a candidate mechanism for 
several cases of bypass transition wherein the 
route to laminar-turbulent transition bypasses 
the well-known paths to turbulence via modal 
instabilities of the underlying laminar basic 
state, for example: 


► large bluntness cones, 

► spherical forebodies, ... 

1 M.V. Morkovin, E. Reshotko, and T. Herbert. "Transition in open flow systems - A reassessment". In: Bull. Am. Phys. 
Soc. 39:1882 (1994). 
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Motivation 


Large bluntness cones 
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Rn=0.15-mm 


Rn = 0.2" 
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Rn=9.53-mm 

Rn = 0.5" 

Rn=12.7-mm 
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Rn = 2.0" 

Rn=50.8-mm 

Small blutness: the transition front moves back relative to sharp cone. 
Large blutness: the transition front moves forward. Why? 


'c.C. Marineau et al. Mach 10 boundary-layer transition experiments on sharp and blunted cones. AIAA Paper 2014-3108. 
2014. 
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Motivation 


Spherical forebodies: the blunt body paradox 


3,4 



Re = 3.0E6/ft Re = 3.9E6/ft Re = 5.0E6/ft 



Figure 14. CEV Global Heating Distribution from Test 6944, a = 28-deg, 8-in. diam. model 


9 Transition observed while the flow is T-S stable, Gortler stable (convex 
curvature) and crossflow is negligible. Why? 


^ Note that images for Re=5.0E6/ft and Re=5.8E6/ft are interchanged. 

^E. Reshotko and A. Tumin. "The blunt body paradox - A case for transient growth". In: Proc. of the IUTAM 
Laminar-Turbulent Symposium V. ed. by H. Fasel and W. Saric. Sedona, AZ, USA, 2000, p. 403. 

4 B.R. Hollis. Blunt-body entry vhicle aerothermodynamics: transition and turbulence on the CEV and MSL configurations. 

AIAA Paper 2010-4984. 2010. 
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Motivation 


Can TG theory explain transition in these cases ? 7 


® Reshotko and Tumin 5 suggested transient growth as a possible explanation 
for transition over blunt nose tips and showed that a transition criterion 
based on linear optimal growth provided successful correlation with the 
PANT database ( M < 6). 

a Recent data 6 corresponds to higher Mach numbers ( M > 6). 


Current ongoing work: 


® supersonic and hypersonic Mach numbers 

9 nonlinear effects (nonlinear streak development, streak instabilities, interaction 
with boundary layer instabilities, ...) 

® spatial approach including non-parallel effects 

9 realistic geometries 

— 

C ' E. Reshotko and A. Tumin. "Role of Transient Growth in Roughness-Induced Transition". In: AIAA J. 42 (2004), 
pp. 766-770. 

® E.C. Marineau et al. Mach 10 boundary-layer transition experiments on sharp and blunted cones. AIAA Paper 2014-3108. 
2014. 

7 P. Paredes et al. Transient growth analysis of compressible boundary layers with parabolized stability equations. (Submitted 
to SciTech 2016). 2016. 
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Zero-pressure-gradient flat plate BL at Mach 3 

® Linear optimal disturbances 

► Verification against prior results using self-similar base flow solution 

► Extension to N-S base flow solution including leading edge shock 
® Non-linear evolution of finite-amplitude linearly optimal disturbances 

® Instability analysis of the perturbed flow 

« Examples of supersonic/hypersonic vehicles with flat plate components: 


® BOEING X-51A (Mach 5) 8 


® NASA X-43A (Mach 10) 9 







Theory 


Transient growth analysis using PSE 


10 



9 The optimal initial disturbance, q 0 , is defined as the inflow condition at Xq 
that leads to maximum energy amplification or gain, G, up to a specified 
position, xi. 

9 Objective function: _/( q) = G = 

o Mack’s energy norm: E(£) = f n q(x) H Mq(x) dQ, with 


M = diag 


f(x) 

jp(x)M 2 


,p(x),p(x),p(x), 


P{x) 

7(7-1 )T(x)M 2 _ ' 


9 Variational formulation using direct and adjoint PSE, 


£(q,q*) = J(q) - (q*,Lq), 


where q* denotes the vector of adjoint disturbance variables. 

9 The PSE approach allows the analysis of a broader class of mean flows 
including those based on NS solvers. 

“j.O. Pralits et al. “Optimal disturbances in three-dimensional boundary-layer flows". In: Ercoftac Bull. 74 (2007), 
pp. 23-31. 
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Theory 


Non-linear development of optimal perturbations 


a The linearly optimal disturbance is used as inflow condition with finite 
amplitude. 

a The nonlinear-perturbation form of the parabolized Navier-Stokes equations 
are used, which is equivalent to the nonlinear plane-marching PSE 11 for a 
single stationary perturbation (A/ = 0, a$ = 0): 


q (x,y,z, t ) 


N 

X ] q n (*,y,z)exp 

n——N 



a An implicit marching technique is used to facilitate the computation of 
nonlinear streaks with large amplitudes, 
a The instability of the resulting modified boundary layer flow is investigated 
by the linear form of the plane-marching PSE 12 . 


11 P. Paredes et al. "The nonlinear PSE-3D concept for transition prediction in flows with a single slowly- varying spatial 
direction". In: Procedia IUTAM 14C (2015), pp. 35—44. 

P. Paredes. “Advances in global instability computations: from incompressible to hypersonic flow". PhD thesis. 
Universidad Politecnica de Madrid, 2014. 
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Results 


Optimal disturbances in flat plate boundary layers 


Verification: Optimal gain vs. spanwise wavenum 


ber 13 1 


a Compressible, self-similar, flat plate boundary layer 

a ReL = 10 s to enable comparison with asymptotic optimal growth results 
based on boundary layer equations 


a M = 1(T 3 


a M = 3, To = 333 K 
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■^P. Andersson, M. Berggren, and D.S. Henningson. "Optimal disturbances and bypass transition in boundary layers". In: 

Phys. Fluids 11 (1999), pp. 134-150. 

^A. Tumin and E. Reshotko. “Optimal Disturbances in Compressible Boundary Layers". In: AIAA J. 41 (2003), 
pp. 2357-2363. 
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Results 


Optimal disturbances in a Mach 3 flat plate boundary layer 


Self-similar vs. full Navier-Stokes solutions 


® Optimal perturbations are commonly computed from the leading edge, 
xo = 0, using the self-similar base flow solution, which ignores both the 
viscous-inviscid interaction near the leading edge and the weak shock wave 
emanating from that region 

® To assess the effects of this simplification, a realistic basic state based on 
the Navier-Stokes (NS) equations is used 


® NS parameters: R„ = 1 /im, 

Re' = 10 6 /m, T 0 = 333 K, T w /T ad = 1 


® Self-similar vs. NS b.l. profiles 

(R = s/xUjv) 





G/Re L 


Results 


Optimal disturbances in a Mach 3 flat plate boundary layer 

Effect of initial optimization position, xq 

» The final optimization position is R\ = yJx\U /v = 1000 with x\ = 


L = 1 m, 


which is used as reference length 

9 Optimal gain using NS base flow 


9 Self-similar vs. NS results 




9 Optimal f3 = 0.25 for inflow location at the leading edge 
9 The overall maximum G / Rei_ is found for xo = 0.25 and = 0.30 
« As expected, the optimal G/ReL results using self-similar and NS base flows agree 
for xo > 0.25, while there is up to 10% deviation for xq closer to the leading edge 
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Results 


Non-linear development of optimal disturbances 


Non-linear evolution of streaks initiated near xq 


0 



9 Ro = \Jx.§U /v — 20 and /3 = 0.25 

« Initial amplitudes defined as Aq — a/£o = A / \J G max with G max — 2398 
9 Streak amplitude based on u: 

As u (x) = [max y z (D(x,y,z)) — min y z (u(x, y, z))] /2 


9 Streak amplitude, As u 


9 Energy gain, G = E / Eq 
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Results 


Non-linear development of optimal disturbances 


Structure of streaks during non-linear evolution 
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Results 


Non-linear development of optimal disturbances 


Non-linear streaks initiated at optimal xq 


9 Ro = \JxqU jv = 500 and /3 = 0.30 


9 Streak amplitude, As u 9 Streak amplitude, As u , for xq — > 0 



R= x fdJJT, R= s /WJ7, 
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Results 


Non-linear development of optimal disturbances 


Structure of streaks initiated at optimal xq during 
non-linear evolution 
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Results 


Instability analysis of the perturbed BL flow 


Instability of streaks initiated at the leading edge 


9 R 0 = sj xqU/v = 20 and (3 = 0.25 

9 Spatial PDE based analysis at the final optimization position, Ri = 1000 
9 Two types of shear layer modes become unstable, the sinuous (S) and the 
varicose (V) modes 


9 Growth rate, —a;, vs. frequency, u 



9 Sinuous, 

4 = 6 , 
u> = 0.115 


9 Varicose, 

4 = 6 , 
u = 0.0958 
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Results 


Instability analysis of the perturbed BL flow 


Amplitude evolution of streak instabilities 


® Ro = s/xqUJu = 20 and /3 = 0.25 
® Plane-marching PSE analysis of the sinuous mode 
® N-factor based on Mack's energy norm 
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Results 


Bypass transition in a Mach 3 flat plate BL 


Natural transition in the unperturbed BL without streak 


® Transition prediction based on Ne — 9 and 10 
® Oblique first mode Ne computed by linear PSE 


® N-factors based on E ® Transition location, Re x = R 2 



R — \JxU jv E 0 = -^Eq 


® N = 10: Re x = 1.0417 x 10 7 , 0 = 0.040, w « 0.012 
® N = 9: Re x = 8.3625 x 10 6 , 0 = 0.045, lu w 0.012 
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Results 


Bypass transition in a Mach 3 flat plate BL 

Bypass transition in a Mach 3 flat plate BL 

® Transition prediction based on Ne = 9 and 10 
a Sinuous mode Ne computed by linear plane-marching PSE 

9 As u vs. £q = 2/ L z Eo 9 Transition location, Re x = R 2 




Eo = £E a 


9 R 0 = 20: P = 0.25, w = 0.075 - 0.100 
9 Rq = 500: p = 0.30, w = 0.075 - 0.100 
« Oblique first mode: j3 = 0.040 — 0.045, lo fa 0.012 




Summary and Conlusions 


® Verification of newly developed linear, adjoint-PSE based optimization code 
for compressible flows with results from literature. 

® The effects of the viscous-inviscid interaction near the leading edge and the 
weak shock wave produce a 10% deviation from results based on self-similar 
base flow when the initial optimization position is located near the leading 
edge. 

® Sinuous mode becomes unstable before the varicose mode for the moderate 
amplitude streaks studied. 

® A non-linear perturbation form of the parabolized Navier-Stokes equations is 
demonstrated for the computation of the downstream development of 
high-amplitude stationary disturbances. 

® The optimal initial optimization position for maximum energy also leads to 
earlier/stronger growth of secondary instability and, hence, earlier onset of 
bypass transition relative to other inflow locations. 

® Bypass transition via finite-amplitude optimal disturbances is demonstrated 
for a zero-pressure-gradient flat plate boundary layer. 
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Extra charts 


Instability of streaks initiated at the optimal xq 

« Ro = \JxqU / v = 500 and /3 = 0.30 

9 Spatial PDE based analysis at the final optimization position, Ri = 1000 
9 Two types of shear layer modes become unstable, the sinuous (S) and the 
varicose (V) modes 


9 Growth rate, —a/, vs. frequency, u> 



9 Sinuous, 
A = 5, 
uj = 0.11 
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Extra charts 


Amplitude evolution of streak instabilities 


9 R 0 = ^xqU/u = 500 and (3 = 0.30 

9 Plane-marching PSE analysis of the sinuous mode along x 
9 N-factor based on Mack's energy norm 
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